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ABSTRACT 


This thesis presents a computer simulation model for 
studying take-off schedules at Kalamata Air Force Base in 
Greece. Six aircraft take-off schedules were examined and 
a comparison of resultS waS based upon factors of 
performance and efficiency/safety. The overall simulation 
model can be easily modified to examine other aircraft take- 


off schedules. 


Lat a 


TABLE OF CONTENTS 


rt INTRODUCTION cee ee ee 


THE PROBLEM sees ee ene 


II. MODEL BEVELORMEN geese) 0c ecm 


Ass DESCRIP TON cece te 


B. MODEL LOGIC Vee 


C. ENPUT AND OUTPUT ep c- eee 


TIll. GPSS APPLICATION eee 


A. BRIEF GPSs DESCR? el 1 Cl peer 


B. MODEL STRUCTURE SINGGES see 


C. GATHERING STATISTICS WITH GPSS 


IV. MODEL VALIDAIAZON “e2a8s oe) eee 


Nite RESULTS—-ANALYSIS We.) are eee 
A. SCHEDULES” ..:) @-e-e ae ee 
B. -RESULTS 3°93 ees) oe eee 


Co SANALYSITS (2 -ee,  eeee 


VI. CONCLUSIONS —= SUMMARY 33... 
As. CONCLUSIONS 2. so voces eons 


B.. “SUMMARY <. 22202 = 5) 0) ne em 


iv 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


49 
54 
60 
62 


73 





I. INTRODUCTION 


THE PROBLEM 
Air Traffic Control has become a safety issue of great 

importance during the last decade because of the many 
near-miss or tragic accidents that have occurred at military 
and civilian airports worldwide. The main reasons for these 
accidents are: 

* Air Traffic Control system failure, 

* Air Traffic Control erroneous procedures, 

ZeePilot error, 

* Weather conditions, 


ea inereased air traffic. 


This issue has even greater significance at military air 
training bases because of the very high ratio of student 
pilots to experienced pilots that are using the air space 
and because of the large volume of aircraft activity in the 
air in specific areas. It is standard military training 
practice for a large number of training aircraft to be 
assigned to the same radio channel for aircraft-to-aircraft 
and for aircraft-to-air traffic controller communications 
and to have a large number of training aircraft following 
approximately the same air pattern with the same air speed 


and altitude. These conditions can increase the probability 


of breakdowns in synchronization and communication between 
aircraft and between the aircraft and the air traffic 
controllers and can therefore increase the probability of 
accidents occurring. 

Kalamata Air Force Base iS an Air Force training base 
located in southern Greece. This Base utilizes two types of 
aircraft for training purposes, T-37's and T-2's. The two 
different types of aircraft used at Kalamata Air Force Base 
ordinarily have different flight schedules, flight 
capabilities, and training missions that would allow them to 
take-off and begin their training flights at different 
times. However because of current operational constraints it 
is common for a group of T-2 and T-37 aircraft to complete 
their missions at approximately the same time so that they 
Simultaneously return to the local Air Traffic Pattern. It 
is during the simultaneous approach of the returning 
aircraft that a critical safety problem arises due to the 
traffic congestion and pilot fatigue. Furthermore, the 
inefficient aircraft schedules interfere with the 
performance of the scheduled activities. A more efficient 
scheduling of aircraft take-offs can remedy the safety 
problems and also can diminish the necessity for aircraft to 
wait in order to get into the mission areas. 

This thesis provides a computer simulation model 
programmed in GPSS for the IBM-PC for analyzing current 


procedures for the efficient management and control of the 


air traffic of Kalamata Air Force Base. The thesis COntcalLns 


the model logic, the GPSS program, model validation and 


model results. 


II. MODEL DEVELOPMENT 


A. DESCRIPTION 
The environment that is modeled is the operation of the 

Kalamata Air Force Training Base in Greece. This operation 
is modeled as a sequential multiserver limited capacity 
queuing system. The server elements of the model consist of 
the base, two runways, an air traffic pattern and eight 
training mission areas. 

The base, runways and air traffic pattern are illustrated 
in Figure 2.1. The air traffic pattern is modeled as a set 
of sequential servers consisting of two entrance points (EPI 
and EP2), the initial point (IP), the low initial point 
(LIP), the break point and the base Key. EP1 is the 
entrance point for the aircraft returning from a western 
mission area and EP2 is the entrance point for aircraft 
returning from an eastern mission area. 


The mission areas are represented in the model as single 


poznts. Specific mission area training activity is not 
modeled. The eight mission areas are the neighboring areas 
around the airport and are illustrated in Figure 2.2. The 


areas west of the airport are numbered one through four and 
the areas east of the airport are numbered five through 


eight. 
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Figure 2.1 Air Traffic Pattern 





Figure 2.2 Mission Areas 


The calling population consists of a basic training 
squadron employing T-37 type aircraft an advanced training 
squadron employing T-2 type aircraft and occasional 
aircraft from other bases. 

B. MODEL LOGIC 
The general flow through the model for a typical aircraft 
is as follows: 

A take-off time is scheduled for the aircraft from one 
of the two squadrons and the aircraft enters the base 
activity queue at that time. If there is no runway or air 
space conflict with other aircraft landing or taking-off, 
the aircraft takes-off and proceeds to the base departure 
point where it requests and is assigned to a mission area to 
Carry out the scheduled training activities. 

If all training mission areas are occupied the aircraft 
is assigned to the mission training area with the smallest 
waiting line and in the case of ties the aircraft is 
assigned to the mission area with the earliest expected 
departure time for the occupying aircraft. All the aircraft 
that wait for a mission area to become available, maintain a 
maximum altitude of 7000 feet at the corresponding areas 
and perform training maneuvers consistent with the altitude 
safety requirements. Once the aircraft enters the mission 
training area it stays in the mission area for a standard 
length of time for the aircraft type and performs the 


scheduled activities. 


At the completion of thE Mmissiom “area trate 
activities the aircraft returns to the base and enters the 
returning air traffic pattern and follows the sequence of 
events as described below and illustrated in Figure 2.1. 

Aircraft returning from the western mission areas enter 
the air traffic pattern at entrance point EP1 and aircraft 
returning from the eastern mission areas enter the air 
traffic pattern at entrance point EP2. The sequence of air 
legs in the base air traffic landing pattern and the 


aircraft capacity of each air leg is as follows. 


Air Leg Aircraft Capacity 
EP1I-EP2 2 
BeZ aie 2 
ELP=LIP 2 
LIP-Break Point 3 
Break Point-Base Key g 
Base Key-Final 2 


If any of the first four air legs is at full capacity 
when an aircraft attempts to enter, the aircraft must orbit 
at the air leg entry point and wait for entry. If the Break 
Point-Base Key leg is at full capacity when an aircraft 
attempts to enter, the aircraft attempting entry returns to 
the LIP-Break Point air leg entry point. If the Base Key- 


Final leg is at full capacity when an aircraft attempts to 


enter, the aircraft attempting entry performs "go around", 
that is, it cancels the landing and attempts a "close 
approach" pattern. 

Aircraft in the "close approach" pattern re-enter the air 
traffic pattern at the Base Key point if no other aircraft 
is in the Break point-Base Key air leg and if no other 
aircraft is waiting for take-off. Otherwise they re-enter 
the air traffic pattern by joining the waiting line at the 
LIP entry point. 

Upon being cleared for landing the aircraft will either 
land and return to the squadron or if flying constraints 
allow the aircraft will perform a touch and go landing. If 
a touch and go landing is made the aircraft will either 
enter the "close approach" pattern or attempt to re-enter 
the air traffic pattern with associated probabilities of 
0.20 and 0.80. If there is a waiting line at the LIP point 
or if the IP-LIP air leg is at full capacity the touch and 
go aircraft will either continue to LIP or re-enter the air 
traffic pattern at the EP1l entry point with associated 
probabilities of 0.70 and 0.30. 

It is also possible that the flow of aircraft in the air 
traffic landing pattern can be interrupted by emergency 
events. Emergency events such as engine failure, oil 
pressure failure, low fuel, or landing gear failure, are 
common occurrences that require the aircraft with the 


emergency to land as soon as possible. When such an 


emergency event occurs, the aircraft with the emergency pre- 
empts all other aircraft in the landing pattern. During the 
emergency the nonemergency aircraft in the landing pattern 
upon reaching the IP, LIP, or Break air leg points leave the 
pattern and return to the entrance point EP1 and join the 
orbiting waiting line. If a nonemergency aircraft is on the 
Base-Final air leg when the emergency Situation arises the 
aircraft continues and makes a full stop landing. 
During the emergency the aircraft present in the 
mission areas do not leave their areas but after completing 
their scheduled activities orbit at a lower altitude until 
the emergency ends to avoid interfering in the local traffic 
pattern. Aircraft returning to the air traffic pattern will 
not enter the pattern but will orbit at the entrance point 
until the emergency ends. 
If an aircraft take off is scheduled during an emergency 
event the take-off is delayed until the emergency event is 


completed. 


ce. INPUT AND OUTPUT 

All input data for the model as well aS suggestions 
for model structure was provided by the instructor 
pilots and from the control tower personnel at the 
Kalamata Air Force BaSe. The general inputs to the model 
consisted of individual aircraft characteristics and 


performance data, alternate take off schedules, mission 


10 


area constraints, time distributions for assignments in the 
mission areas, and the time distances between the 
reference points for the air legs in the air traffic 
pattern. 


The model outputs contain information pertaining to the 


aircraft of each squadron, the air traffic pattern and 
mission area utilization. The following output is 
available. 


fMm—eriticetrartt total flayqht time distribution. 
2. Number of aircraft by squadron in the east and 


the west mission areas. 


3. Number of aircraft take-offs. 
4. Number of entries in each mission area. 
5. Number of entries in each leg of the air 


traffic pattern. 
6. Maximum number of aircraft waiting in each 
queue. 
7. Average wait time for each queue. 
8. Waiting time distributions for the following 
queues: | 
ae beet larric Controller. 
b. Entry points, EP1 and EP2. 
Co imtial point, IP. 


ae LOw initial point, BLiP. 


ite 


III. GPSS APPLICATION 


A. BRIEF GPSS DESCRIPTION 
GPSS is the General Purpose Simulation System language 
developed by IBM for modeling and simulating queuing 
systems. GPSS was used to model the Kalamata Air Force Base 
air traffic operations. The GPSS program is included in 
Appendix A. 
GPSS uses the process interaction approach for modeling 
in which the model entities are either temporary or 
permanent. The temporary entities are called 
transactions and the permanent entities are called 
facilities and storages. The transactions represent the 
calling population and the facilities and storages 
represent the service centers. Transactions interact with 
other transactions and with the facilities and 
storages. In the Kalamata Air Force Base model the 
calling populations of aircraft are represented by 
transactions and the mission areas and the air traffic 
landing pattern segments are represented by facilities 
and storages. 

The modeling and programming approach in GPSS is to 
define a set of programming statements called blocks 
that represent the entrance and flow of the 


transactions into the queuing system composed of the 
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facilities and storages. There can be many 
transactions Simultaneously moving through the blocks. 
At any point in time each transaction is positioned at a 
block and most blocks can hold many transactions 
Simultaneously. The transfer of a transaction from one 
block to another occurs instantaneously at a specific time 
or when some change of system condition occurs. Time in 
the GPSS model is managed by the next event sequence 
with the simulation clock changing at nonuniform discrete 
time points when the state of the system changes. 
Transactions continue to move through the system until they 
either encounter a waiting line or service time delay. 

In GPSS simulated clock time is an integer value 
whose scale value is chosen by the programmer. The unit 
of time is not specifically stated but is implied by 
providing all times in the same units. In the Kalamata Air 


Force Base model the unit of time used is the second. 


B. MODEL STRUCTURE IN GPSS 

A brief description of some of the important 
programming blocks and storage areas used in the GPSS 
Kalamata Air Force Base model are contained in this 
section. 

GENERATE and TERMINATE blocks: Transactions are created 
and enter the system at one or more GENERATE blocks and 


are removed from the Simulation at TERMINATE blocks. 


1S 


The time and frequency with which transactions enter the 
system are controlled by the GENERATE block. In the 
Kalamata Air Force Base model GENERATE blocks are used for 
the entry of training aircraft from the squadrons’ for 
take off assignments and for the entry of occasional 
aircraft from other bases into the Kalamata air trateite 
landing pattern. 

ADVANCE block: The ADVANCE block will hold 
transactions for a specified or computed number of time 
units. The purpose of the ADVANCE biock is to _ hold 
the transactions in service. In the Kalamata Air Force 
Base model ADVANCE blocks are used to simulate the time 
delays associated with take off delays, training mission 
areas and transit from point to point in the air Wee 
segments of the air traffic landing pattern. 

TEST block: The TEST block is used to manage or transfer 
transactions based upon the test conditions. In the 
Kalamata Air Force Base model the TEST block is used to 
prevent aircraft from entering the system after the daily 
training period and to assure that the aircraft in the 
system are correctly processed in order to complete all 


landings after the daily training period ends. 


GATE block: The GATE block is used as a gate to 
interrupt the flow of transactions depending upon 
conditions that set the gate to "open" or "closed". In the 


Kalamata Air Force Base model GATE blocks are used to 
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prevent aircraft from continuing in the air traffic landing 
pattern or the take-off queue during an emergency event. 

SELECT block: The SELECT block is used to direct the 
flow of transactions. In the Kalamata Air Force Base 
model SELECT blocks are used to assign the aircraft to the 
training mission areas after take-off. 

JOIN, REMOVE, COUNT and MARK blocks: The JOIN, REMOVE, 
COUNT and MARK blocks are used to collect, remove, 
count and identify transactions in the queue. In the 
Kalamata Air Force Base model if the training mission areas 
are occupied the JOIN, REMOVE COUNT and MARK blocks 
are used to determine current aircraft assignments based 
upon the shortest waiting lines for the areas. 

FACILITIES and STORAGE areas: GPSS FACILITIES and 
STORAGE areaS are used to collect and hold transactions 
for time delays that can be associated with service or 
performance of the transactions. A GPSS FACILITY can 
hold one transaction. A GPSS STORAGE area can hold more 
then one transaction. In the Kalamata Air Force Base 
model eight FACILITIES model the eight mission areas, one 
FACILITY models emergency aircraft, One FACILITY models 
the air tCrartic¢e controller, four facilities model 
aircraft synchronization, and six STORAGE AREAS have been 
used to model the air leg segments of the air traffic 


landing pattern as described in the following page. 


ILS: 


FACILITY 


IO MIL 
LOZ 
Tes 
104 
IC SIs) 
106 
tO? 
108 
DANGER 
CNTR 
DUMY1 
DUMY 2 
DUMY3 


DUMY4 


STORAGE AREA 


CAP1: 
CAPZ: 
CAP3: 
CAP4: 
CAP5: 


CAP6: 
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MODEL 


Mission area l. 
Mission area 2 
Mission area 3. 
Mission area 4 
Mission area 5. 

Mission area 6. 

Mission area 7. 

Mission area 8. 

Emergency aircraft. 
Traffic controller. 
Aircraft synchronization. 
Aircraft synchronization. 
Aircraft synchronization. 


Aircraft synchronization. 


Air Leg Segment 


EPRAEPZ 

EP2-IP 

IP-LIP 

LIP-BREAK POINT 
BREAK POINT-BASE KEY 


BASE-FINAL 


e. GATHERING STATISTICS WITH GPSS 

GPSS automatically records data and collects queue 
statistics for transactions that pass through a storage 
area. In addition to the previously described GPSS 
programming blocks that manage the flow of transactions 
there are several block types that are specifically designed 
to gather statistics on transactions. These blocks and 
their application are described in this section. 

QUEUE and DEPART statistic blocks: QUEUE and DEPART 
blocks are used to identify specific data collection 
points in the queue. 

Data on transactions that move through the queue and 
enter and leave associated QUEUE and DEPART blocks is 
collected as separate sets of queue statistics. In the 
Kalamata Air Force Base model queue statistics were 
accumulated over the entire period of the simulation, in 
each queue of the air traffic pattern and in each queue of 
the mission areas. The accumulated statistics for each of 


the model queues are identified in the following list. 


Queue Statistics 

Maximum number of aircraft. 

Average number of aircraft 

Standard deviation of the number of aircraft. 
Average waiting time. 


Number of entries. 
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Number of entries that did not wait for entry. 


Average wait time to enter the queue. 


TABULATE, TABLE and QTABLE statistic blocks: The 
TABULATE, TABLE and QTABLE blocks are used to collect 
data for frequency and cumulative frequency tables. 
The TABLE and QTABLE blocks define the transaction 
characteristics that are to be counted and the range of 
the frequency tables. Data on characteristics of the 
transactions that pass through a TABULATE block are 
automatically collected. The TABULATE block may be used 
anywhere in the GPSS program and can collect data on all 
of the transactions computed characteristics. The 
TABULATE, TABLE and QTABLE blocks are used in the Kalamata 
Air Force Base model to tabulate the total flight time 
distribution for the two types of aircraft and waiting time 


distributions for all of the model queues. 
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IV. MODEL VALIDATION 


The Kalamata Air Force Base model was validated by 
comparing model output to historical data for a specific 
take off schedule for the squadrons of aircraft. Historical 
Gata was made available for each aircraft for a three day 
period. The historical data is included in Appendix C. 

For each aircraft the historical data consisted of the 
interarrival times for take off, the time spent waiting in 
the air traffic controller queue, the time spent waiting to 
enter area 3 and the total flight time. The historical data 
waS accumulated and averaged for the three days of base 
activity. Using the historical interarrival rates as model 
inputs the model generated waiting times and total flight 
times were compared to the historical data. 

The historical interarrival rate appeared to be 
nonhomogeneous on a daily basis. Therefore the operations 
day was broken into time periods for which the interarrival 


rates were homogeneous. These time periods were: 


1) 07:30 - 08:40 5) 11:40 - 14:00 
2) 08:40 - 09:10 6) 14:00 - 14:30 
3) 09:10 - 11:10 7) 14:30 - 16:20 
4) 11:10 - 11:40 8) 16:20 - 17:30 
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For each time period the interarrival times appeared to 
have an exponential distribution. This hypothesis was 
tested using the Kolmogorov-Smirnov test and was not 
rejected. The hypothesis test results are included in 
appendix D. With this base information the model was run 
for a period of eight days using the hypothesized 
exponential distributions for the take off interarrival 
times. The validation comparison tests follow. 

Figure 4.1 displays a plot of the comparison of the 
actual wait time with model generated wait time for entry to 
area 3, and Figure 4.la displays the regression line for the 
regression model: 

Simulated Data=A+B*(Actual Data). 

The slope of the regression line is .97 with a standard 
Geviation of .027. The hypothesis B=1l was tested with the 
t-test and was not rejected. The probability level for the 
t-statistic was .16. 

Figure 4.2 displays a plot for the comparison of actual 
wait time with model generated wait time in the air traffic 
controller queue and Figure 4.2a displays the regression 
line for the regression model: 

Simulated Data=A+B* (Actual Data). 

The slope of the regression line is .98 with a standard 
deviation of .043. The hypothesis B=1l was tested with the 
t-test and was not rejected. The probability level for the 


t-statistic was .43. 
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Figure 4.3 displays a plot of the comparison of actual 
total flight time with model generated total flight time 
and Figure 4.3a displays the regression line for the 
regression model: 

Simulated Data=A+B*(Actual Data). 

The slope of the regression line is .96 with a standard 
deviation of .0O19. The hypothesis B-1 was tested with the 
t-test and was not rejected. The probability level for the 
t-statistic was .09. 

These results provide a validation of the model and show 
that the model data compares favorably with the historical 


data for the specified take off schedules. 
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Figure 4.1 Comparison of Actual Wait Time 


with Simulated Wait Time in Area 3 
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Figure 4.la Regression of Simulated Wait Time 


on Actual Wait Time in Area 3 
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Figure 4.3 Comparison of Actual Flight Time 


with Simulated Flight Time 
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V. RESULTS~ANALYSIS 

This section contains a description of four basic and 
two additional take-off schedules for the squadron of 
training aircraft for Kalamata Air Force Base that were run 
in the model and a comparison analysis of the resulting air 
traffic controller congestion in order to determine a 
reasonable efficient take-off schedule. In order to @Eua 
each take-off schedule in the model, the GPSS program had to 
be modified. These program modifications are included in 


Appence 06. 


A. SCHEDULES 

Schedule 1: This schedule consists of a forty minute 
take-off cycle. The cycle consists of two consecutive 
fifteen minute take-off periods followed by ten minutes of 
no take-off activity. The T-37 aircraft squadron assigns 
five aircraft for take-off in the first fifteen minute 
period. The T-2 aircraft squadron assigns five aircraft for 
take-off in the second fifteen minute period. This cycle is 
repeated until all the training aircraft are scheduled for 
take-off. Aircraft that can not take-off in their assigned 
period are recycled for later take-off. Aircraft taking-off 
have priority use of the runway over normally landing 
aircraft. This schedule was proposed for analysis by the 


two squadrons at Kalamata Air Force Base. 
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Schedule 1A: This schedule is derived from and is 
identical to Schedule 1 except that the number of aircraft 
scheduled for take-off in each fifteen minute period is four 
instead of five. This schedule was suggested as a result of 
reviewing the model output data for Schedule 1. It was 
thought that this change would decrease the number of 
training flights while also decreasing the air traffic 
congestion and waiting times. 

Schedule 1B: This schedule is derived from and is 
identical to Schedule 1 with the take-off cycle period 
extended to sixty minutes. This schedule was also suggested 
aS a result of reviewing the model output data for Schedule 
i. It was thought that this change would also decrease the 
number of training flights while decreasing the air traffic 
congestion and waiting times. 

Schedule 2: This schedule consists of an eight-minute 
take-off cycle. The cycle consists of two consecutive 
four-minute take-off periods. The T-37 aircraft squadron 
assigns one aircraft for eco Oe in the first four-minute 
period. The T-2 aircraft squadron assigns one aircraft for 
take-off in the second four-minute period. This cycle is 
repeated until all ene training aircraft are scheduled for 
take-off. Aircraft that can not take-off in their assigned 
four-minute period are recycled for later take-off. MThis 
schedule was also suggested as a result of reviewing the 


model output data for Schedule 1. It was thought that this 


ARS, 


model output data for Schedule 1. It was thought that this 
change would increase the number of training flights while 
decreasing the air traffic congestion and waiting times. 
Schedule 2A: This schedule is derived from Schedule 2 and 
is identical to Schedule 2 except that the take-off cycle is 
extended to ten minutes with two five-minute take- off 
periods. This schedule was suggested as a result of 
reviewing the model output data for Schedule 2. 

Schedule 3: This schedule is a reproduction of the 
schedule currently in use at Kalamata Air Force BaSe. The 
schedule contains no structure and take-off times are 
scheduled at random. The distributions o£ CUpmere 
interarrival times for take-off times were analyzed and 
found to be exponential for different periods during the 
day. This analysis is contained in Appendix C and was also 


used for model validation. 


B. RESULTS 

The model output data for each schedule are contained in 
Appendix E. The measures of effectiveness used in comparing 
schedules were performance and efficiency/safety. 
Performance iS measured by the average number of training 
aircraft scheduled. Efficiency is measured by the average 
waiting times in the mission area queues. Safety is 
directly related to efficiency in that the smaller waiting 


times mean less hazardous flying conditions. The empirical 
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Gistributions for both of these measures for each schedule 
are contained in Appendix E. 

The summary results for Schedule 1 are of interest because 
this schedule was suggested by the training squadrons at 
Kalamata Air Force Base. This schedule, by comparison with 
the other schedule results, does not provide the highest 
values for performance and efficiency/safety. It was 
because of these results that the other schedules were 
derived from Schedule 1 by sensitivity analysis. 

The summary results for Schedule 3 are of interest 
because this schedule is a reproduction of the schedule 
currently in use at Kalamata Air Force Base. This schedule 
by comparison with the other schedule results, also does not 
provide the best values for performance and 
efficiency/safety. 

It was found that by applying Schedule 1A and 
comparing the results to Schedule 1 that performance 
decreased but that efficiency/safety improved drastically. 
Figure 5.1 compares the cumulative mission area waiting time 
distributions and Figure 5.2 compares the cumulative LIP 
point waiting time distributions for Schedules 1 and 1A. 

Schedule 2 results were more efficient than Schedule 1 
results. However, by applying Schedule 2A mission area 
waiting times decreased dramatically and provided the most 
efficiency/safety measure for all schedules examined. 


Figure 5.3 compares the cumulative mission area waiting time 
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distributions and figure 5.4 compares the cumulative LIP 
point waiting time distributions for schedules 2 and 2A. 
Figure 5.5 compares the cumulative mission area waiting time 
distributions for schedules 1A, 2A and 3. It is obvious 
that Schedule 2A is preferable to Schedule l. 

Schedule 3 results were more efficient than Schedule 1 
results. The performance measure for Schedule 3 however is 
the lowest of all the schedules examined. Figure 5.6 
compares the cumulative LIP point waiting time distributions 


for schedules 1A, 1B, 2A and 3. 


C. ANALYSIS 

From the above summary results the preferred schedules 
appear to be Schedules 1A and 2A. An analysis of vardemee 
was performed, using the function "ANOVA" from the OA3660 
APL WORKSPACE, to test the hypothesis that the mean mission 
area waiting time differences for schedules 1, 1A, 2, 2A are 
not significant. The analysis of variance results are 
contained in Table 5.1. These results show that the null 
hypothesis of no significant differences between mean 
mission area waiting times is rejected at significance 
levels greater than .995. 
The Sum of Squares from the previous analysis was broken 
into three components in order to test for individual 
effects rather than just a schedule effect using individual 


degrees of freedom. The results of this analysis are 
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contained in the tables 5.2 and 5.2a. These results show 
that the hypothesis of no significant difference between the 
compared mean mission area waiting times for the selected 
schedules is rejected for each comparison at significance 


levels greater than .975. 


TABLE 5.1 


ANALYSIS OF VARIANCE RESULTS 


ANOVA TABLE 


SOURCE 


SCHEDULE 252864.74 Va Ee 7, 


ERROR Pei. OS CIT D2 7. 
TOTAL 7 244045.78 

R-SQUARE = 0.954 

OVERALL MEAN = 200.68 

TREATMENT EFFECTS -87.56266.68%**.50 14.39 
oe = 2 : Onasor -17 2.92 


BS wo = 2400 2 =O75 S59 “17.92 
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TABLE 5.2 


ANOVA WITH INDIVIDUAL DEGREES OF FREEDOM 


SCHEDULE 1 SCHEDULE 1A SCHEDULE 2A SCHEDULE 2 


44.62 181.62 49.138 443.34 7202.0 .7 197.15 232593 
1 1 al 1 1 li 1 1 


716.45 64162.575 
-708.48 125485.970 
415.77 4A32Z16.1732 
=o 200 9384.500 
48.04 EES 31920 
Oe) 0.378 


1284.50 642.250 
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TABLE 5.2a 


ANALYSIS OF VARIANCE RESULTS 


USING INDIVIDUAL DEGREES OF FREEDOM 


(1,1A)vs(2, 2A) 64162.575 64162.575 


1 vs 1A 125485 .970 125485 .970 
2vs 2A BOZO. 4 7 3 43216.173 


RESIDUAL 11181.048 2795.262 





35 


—. 


— 


tot a wet © et 


ot oF oe ee 


SUMEUULES 1 AND 1A vs WHITING TIME 





es aay 
e a SH 


VRITING TIME 


Figure 5.1 Cumulative Waiting Time Distribution 


in the Mission Areas 


Schedules 1 and 1A 
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Figure 5.2 Cumulative LIP Waiting Times 


Schedule 1 and 1A 
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Figure 5.3 Cumulative Waiting Time Distribution 
in the Mission Areas 


Schedules 2 and 2A 
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Figure 5.4 Cumulative LIP Waiting Times 


Schedule 2 and 2A 
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Figure 5.5 Cumulative Waiting Time Distribution 
in the Mission Areas 


Schedules 1A 2A and 3 
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Figure 5.6 Cumulative LIP Waiting Times 
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VI. CONCLUSIONS - SUMMARY 


A. CONCLUSIONS 

If an added emphasis 18s to be placed upon the 
efficiency/safety factor in scheduling the aircraft then 
Schedule 2A 1S the preferred schedule. If Schedule 2A 
conflicts with local base operations due to other operating 
constraints then Schedule 1A is the next best schedule. The 
difference in mean mission area waiting times for these two 
schedules is made more significant when it is also 
realized that landing aircraft have runway priority in 
schedule 2A and do not in Schedule 1A. Figure 6.1 compares 
the cumulative mission area waiting time distributions and 
figure 6.2 compares the cumulative LIP point waiting time 
distributions for schedules 1A and 2A. | 

If added emphasis is to be placed upon the performance 
factor in scheduling the aircraft then Schedule 2 is the 
preferred schedule. However, the increase of twenty four 
scheduled aircraft is at the expense of more than a ten fold 
increase in mean mission area waiting times that contributes 
to air traffic congestion and pilot fatigue and should be 
avoided if possible. Again, if Schedule 2 conflicts with 
local base operations due to other operating constraints 


then Schedule 1 is the next best schedule. It is 
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recommended that the efficiency/safety factor be the 
deciding factor in selecting schedules. 

Table 6.1 contains a summary comparison between the 
different schedules. In this table the first column 
contains the average number of flights, the second column 
contains the percentage of the aircraft that waited in the 
mission area queues, the third column contains the 
percentage of the aircraft that waited more than 180 seconds 
in the mission area queues, the fourth column contains the 
percentage of the aircraft that waited more than 30 seconds 
at LIP and the fifth column contains the average conditional 


waiting time, in seconds, in the mission area queues. 


TABLE. 6. 


SUMMARY COMPARISON BETWEEN SCHEDULES 


# of A/C |% of A/c!|% of A/C/% of A/C|Avg.Conditional 
waiting time 
in seconds 


Mm. areas 





Figure 6.3 compares the cumulative mission area waiting 
time distributions and figure 6.4 compares the cumulative 


LIP point waiting time distributions for schedules 1 and 2. 
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Figure 6.1 Cumulative Waiting Time Distribution 
in the Mission Areas 


Schedules 1A and 2A 
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Figure 6.2 Cumulative LIP Waiting Times 


Schedules 1A and 2A 
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Figure 6.3 Cumulative Waiting Time Distribution 
in the Mission Areas 


Schedules 1 and 2 
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B. SUMMARY 

This thesis presents an IBM-PC GPSS model for studying 
aircraft take-off schedules at Kalamata Air Force Base in 
Greece. The model is specific to Kalamata Air Force Base 
because of the structure and sequence of the queues used in 
modeling the air traffic control points and the training 
mission areas. The GPSS program can be transferred to any 
IBM-PC compatible computer that will run GPSS. The model 
can therefore be used at the Kalamata Air Force Base to 
continue to examine other aircraft take-off schedules and 
to help reducing fuel consumption. 

Six aircraft take-off schedules were examined and a 
comparison of results was based upon factors of performance 
and efficiency/safety. The take-off schedule currently used 
at Kalamata Air Force Base is shown to be a poor performer 
with a low efficiency/safety factor. The take-off schedule 
proposed for use by the training Squadron personnel at 
Kalamata Air Force Base is shown to have the highest 
performance factor and the worst efficiency/safety factor. 

The specific schedule of ten minutes of a ten minute take- 
off cycle consisting of two consecutive five minute take-off 
periods from each squadron with runway priority given to 
landing aircraft is shown to be the best schedule based upon 
the recommended emphasis placed upon efficiency/safety 


rather than on numerical performance. 
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R260 4 MeAuie? TG BR le ce 
22/0 
2286 TOTALZ? TABLE MP7 ,3400,240,19 *TOTAL TIME DISTRIBUTIGN 


2270 TOTALT2 TABLE MF6,5500,240,18 ;TOTAL TIME DISTRIBUTION 


2300 CNTR Q@TABLE CNTR,O, 60,40 -TIME DISTRIBUTION IN CNTR QUELIE 
2310 CAP1 TABLE CAF1,0,30,40 -TIME DISTRIBUTION IN CAF1 QUEUES 
2320 CAP2 OQTABLE CAP2,0,30,40 sTIME DISTRIBUTION IN CAP2 QUEUE 
2330 CAPZ  QTABLE CAPS ,0,30,40 sTIME DISTRIBUTION IN CAPS QUEUE 
2340 CAP4 QTABLE CAP4,0,30,40 :TIME DISTRIBUTION IN CAP4 QUEL 

2250 AREAL QTABLE 101,0,120,22 -TIME DISTRIBUTION IN AREA Q10: 
2360 AREAZ QTABLE 102-0 120eec sTIME DISTRIBUTION IN AREA Gié2 
2370 AREAS GTABLE 103,0,120,32 -TIME DISTRIBUTION IN AREA Gi0z 
2590 ARTAS OTABLE P04 0. 120,52 “TIME DISTRIBUTION IN AREA 2104 


Mee Gnts= 286 OTASLE Use ol eae Suis DIShe THs tOee in GAES 
meee FRESE CTASLE 106,0,120,32 ;TIME DISTRIBUTION IN ARE 
Pee Shaay | STAELE ee ens! sTine DISTRIBUTION In GRE 
Pee AREAS GTAELE FOG Ws 20.52 PRIME DISTRIBUTION. IN eRe 


2440 FIRST BVARTABLE EO tere Neste E 1OStP2°E°104) 
24609 * ES ae TBE Sle GeneA.& ES 


2470 CAPI STORASE 
2480 CAPZ TORAGE 
2490 CAPS STORAGE 
20900 CAFS STORAGE 
Bee CAPS ShORes= 
2920 CAFS STORASE 


KITA RIPE 


2940 *% ee eee oor ENT 


feee ONT GERSRAIE 29600 sTIMc ARRIVES AT 20606 
20 tesa © N#DOWN NEF INISH sWAIT UNTIL ALL THE A/C LAND 
a7 TERMINATE 1 fonuk GFF ThE RUN 


ae 


ay 


ri 


for the different 


This appendix 


AFFENDIX B 


contains the GPSS proecram meam tea ule. 


Gake--6f i esehned wiecs 


SCHEDULE ©] 


LQ HHH K KEK HKHE HHH HEHEHE HEHEHE KEKE HEHE KEE KEK KEKE EEE KEKE KEKE KEKE KEKKEKHKEK KEKE KHER KEK KEKRKEKRKRHEKEE 
20 * * 
Sa AIR-TRAFFIC CONTROL SIMULATION # 
GO a a ae ea a a es a a a sm ¥ 
oO + + 
EQ FHKE KKEEHKE KEKE HEE KEKE KEEFE FREER EKEEREKEKEKEEKEE EERE KEKE KEEKRERKEEKERKEKEKRKEKEKKEKEKRE KEKE 
7O * 

BO * 

ae EXPONENTIAL DISTRIBUTION 

100 + 

P10 EXE FUNCTION REG, C24 

O,0/.1,-104/.2,.222/,3,. 350 08 SOF! 569/665. 715/.7 1. 2/. 75, Ae 

ee hee 4 2 B.5/,.88,2. 22/.9,2.3/.92 Panes <7 aneet. a Bet Ts sto. ome 

a eee yA 94 5 aie 9) r 99,4. &/, D5 das Say ae 998.6. fy Ae 999 ,7/. 99982 

Oe 

Lee 

140 MEAN FUNCTION fies 


3600, 389.11/5400, 221, 60/1s200, 450,62 


Lo000, 242.6357 25400, 557.927 2o2), | 


24.58 


fi 


21800 , 494. 69/36690, 221.29 


150 * 
160 * 
170 GENERATE ve sCREATE A SINGLE TRANSACTION 
180 SWITCH! LOGIC § { -SET LOGIC SWITCH 1 
190 ADVANCE =00 -T-37 TAKE OFF 
POC LOGIC R { =RESET LOGIC SWICH 1 
710 ADVANCE 2100 >NONE T-37 TAKES OFF 
220 TRANSFER SSWITCH1 
230 GENERATE pee cdn 
a6 6 Menrete (Gene 2 *SET LOGIC SWICH 2 
250) VANCE 200 -T-2 TAKE OFF 
re LOGIC R A *RESET LOGIC SWICH 2 
270 ADVANCE 2100 -NONE T-2 TAKES OFF 
280 TRANSFES SWITCH? 
296 GENERATE 9000 ,FNFEXP :CREATE EMERGENCY EVENTS 
300 SEIZ DANG 
210 ADVANCE 400,180 sEMERGENCY HOLDS 
320) RELEASE DANS 
33 TERMINATE sEMERGENCY TERMINATES 
240 * 
350 GENERATE 9100, FNSEXP 2A/C FROM OTHER BASES 
360 ASSIGN 1,2 ?ASSIGN TO A PARAMETER 
270 MARK 4 =MARK THE TIME 
380 TRANSFER EPI 

% 
409 * CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 

5d 
420 GENERATE 60,10 sCREATE T-37'S 
430 TEST L N¢GNRT,1,PISTA = THE SIMULATION STOPS AFTER 
440 ee be 1,PISTA 
450 ASSIGN 1,1 sASSIGN TO A PARAMETER 
460 TRANSFER DOWN 
A7G + 
42) GENERATE 69.10 reese yaa°e 


100 
110 


: ay . 


oh 


140 
360 
150 
z19 
150 * 
140 
170 
180 
190 
200 
oN ml) “8 
320 
230 
ae 

2@ 


KEEEEFKEEE ERE EEE KEKRFEER EEKEEEFEKEE EERE EEEEEREKKREKRKEKEREKRKEEE EERE EER HER EREEASE 


* + 


¥ AIR-TRAFFIC CONTROL SIMULATION + 
Fr ae ee ee ce ste Sate cate ce cate cate etn ste Stm Stn cate SED St” GED SEED stn nen stn em . 
% Pe 


KEKE EHEKEKEKEKEHEHEKKEKEHEKEKEKEEEHEHEKEEEHEKREEKRKEEEE KKK KKREKRKEKRKEHEKRKEEHEEERKEKHEHRKHETE 


# EXPONENTIAL DISTRIBUTION 
+ 


EXP FUNCTION RN1,C24 
ler enema? 2 fen) 4. HOR / O69 76s 915/.7,5.2/.75, 1.36 
1. a eames .86.2.127.9,2.3/.92,2. ey 94,2.81/.95.2.99/. 54° i 


5 /.98.%.9/.99,4.6/.995,5.3/. 252° &.2/.999,7/.9998,8 
or 

+ 

MEAN = FUNCTION C1,De 


0,380. 11/5400, 221 .60/13200, 488. ree 
00.242, 63/23400, Bre B3/25200,124, 


60,494. 60/3660: 221.3 
Sg 
GENERATE wed SCREATE A SINGLE TRANSAZTION 
SWITCH1 LOGIC & i SSET LOGIC SWITCH 1 
ADVANCE SOC. 21-37 TAKE OFF 
Pere 1 3 RESET LOGIC SWICK 1 
ADVANCE 2100 NONE T-27 TAKES OFF 
TRANSFER SSWITCH! 
GENERATE 13700, 1 
BbTcu? LOGIC S - ;SET LOGIC SWICH 2 
ADVANCE 300 2 TAKE OFF 
POSneUR p i RESET LOGIC SWICH 2 
ADVANTE 7100 “NONE T-2 TAKES OFF 
TRANSFER SWITCH? 
GENERATE 8000 ENEEXF :CREATE EMERGENCY EVENTS 
SEneLe DANGE 
ADVANCE re 186 sEMERGENCY HOLDS 
RELEASE DANCE 
TERMINATE sEMERGENCY TERMINATES 
* 
GENERATE 9100, FNSEXP ;A/C FROM OTHER BASES 
. ASSIGN ae -ASSIGN TO A PARAMETER 
MARK: 4 MARK THE TIME 
TRANSFER EFI 
& 
x CREATE AIRCRAFTS FROM TWO DIFFERENT S@QUADROMS ON BASE 
¥- 
GENERATE 75,10 sCREATE T-37'S 
TEST C NSGNRT,1,PISTA ?THE SIMULATION STOPS AFTER 
GATE LS 1,FISTA 
ASSIGN lel -ASSIGN TO A PARAMETER 
TRANSFER _, DOWN 
% 
GENERATE 75,10 sCREATE T-2'S 
ToS NSGNRT,1,PISTA 
GATE LS 2, PISTA 
ASSIGN ‘0 sASSIGN TO A PARAMETER 
DOWN QUEUE CRTR -WAIT QUT OF THE RUNWAY 
GATE FV DANSE ;DONT MOVE IF EMERGENCY HOLDS 
a GATE SE CAFS =CHECK FOR A\C ON BASE LEG 


on 
cn 


i 
129 
120 
1 


] 


* 
rn 


STE ee 


FEE EEEEREFEREREEFEEREEEEEEEREFREFEEEEEREREREEEEREREEEREREREREEE ERS EPH E 


> * 


¥ 
% 
% 


% 
¥ 
+ 
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AIR-TRAFFIC CONTROL SIMULATION 


KEKEKKEKKEHEKEKEKEEKRKEEREEKREREREREEKEKKEEERKEEKRHEREEEKRHKEKEEEEKRERERKREREREE ER 


% 
¥ 


\. MEAN 
3400. 386.11/5400 
15000. 242. 63/23400, 557 


inne ee 
s 1.6/,. 84 i 
yea See 


FUNCTION 
82/.88, 
.9/.95, 


FUNCTION 


“ 4 
FLIM 


31805 ,494.60/36600, 2 


150 
160 
179 
180 
190 
200 
21¢ 
220 


2350 


EXPONENTIAL DISTRIBUTION 


RN1,C24 
55/24. 509 / no, |G anoueono/ 7 

17 DSc o/ a9 2 2552) 94-2 81) one 

9955, 3/.998, 6. 2/.999,7/.9998,8 


5 1.38 
99/.96,3.2 


° 
= 


Beevers 


ry van 


z 
is 


4.6 


C1,D8 
260/13200, 488.62 

.83/25200,124.58 
Tee 


eae 
Giz 


* 

a 
GENERATE Pie! ;CREATE A SINGLE TRANSACTION 

SWITCH1 LOGIC § i ?SET LOGIC SWITCH 1 
ADVANCE 200 *T=37 TAKE OFF 
LOGIC FR 1 *RESET LOGIC SWICH 1 
ADVANCE 3300 -NONE T-37 TAKES OFF 
TRANSFER ,SWITCH1 
GENERATE occ 

SWITCH? LOGIC S * -SET LOGIC SWICH 2 
ADVANCE 200 -T-2 TAKE OFF 
Lec 2 sRESET LOGIC SWICH 2 
ADVANCE 3300 -NONE T-2 TAKES OFF 
TRANSFER SWITCH? 
GENERATE 9000, FNSEXF sCREATE EMERGENCY EVENTS 
SE ii nle DANG 
ADVANCE 600 , 180 s EMERGENCY HOLDS 
RELEASE DANSE 
TERMINATE sEMERGENCY TERMINATES 

% 
GENERATE 9100, FNSEXF :A/C FROM OTHER BASES 
ASSIGN &D -ASSIGN TC A PARAMETER 
MARK 4 7MARK THE TIME 
TRANSFER EFI 

% 

* CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 

¥ 
GENERATE 60,10 sCREATE T-37'S 
TEST L N¢GNRT,1,PISTA >THE SIMULATION STOPS AFTER 
GATE LS 1,PISTA 
ASSIGN el sASSIGN TO A PARAMETER 
TRANSFER , DOWN 

* 
GENERATE 60,10 sCREATE T-2'S 
TEST L NSGNRT,1,PISTA 
GATE LS 2,PISTA 
ASSIGN 150 sASSIGN TO A PARAMETER 

DS DOWN QUEUE CNTR :WAIT OUT OF THE RUNWAY 

GATE FV DANS= *DONT MOVE IF EMERGENCY HOLDS 

* GATE SE CAFE :CHECK FOR A\C ON BASE LEE 


> 
30 +¥ AIR-TRAFFIC CONTROL SIMULATION 7” 
Be a ae ee ee a a a a a ee ee me + 
50 * ¥ 
EQ HH KK HH HHH EHH EKKEHEK EE KEEKHEKHEE KE KEKE EHEH EERE EHEKEE KKK KKEKEKKKKKE KE KKH EEES 


1 EXPONENTIAL DISTRIBUTION 

Ls 

110 EXP FUNCTION RN1,C24 

1 04/.2,.222/,35,.255/.45.509/.5,.69/.6,.915/.7,1.2/.75 
Seme72G4,1.85/.80,2.127.9,2.5/.92,2.52/.94,2.81/,95,2.99/.96 
697 5 545/. 983. 9/.99,4.6/ 9955.37. 998, 6. 2/.999,7/. 9998 8 

120 * 

1=0 + 

140 MEAN = FUNCTION C1,De 

3600, 380. 11/5400, 221. 60/12200, 488. 62 

15000, 242, 63/23400, 557. 82/ 25200, 124. 58 

%1800,494. 60/24600, 271.20 


150 * 
160 * 
290 GENERATE 9000 , FNSEXF sCREATE EMERGENZY EVEN7S 
200 Sei7E DANGER 
310 ADVANCE 600, 180 sEMERGENCY HOLDS 
320 RELEAS DANGER 
330 TERMINATE sEMERGENCY TERMINATES 
SH) ¥ 
350 GENERATE 9100, FNSEXP sA/C FROM OTHER BASES 
360 ASSIGN ee -ASSIGN TO & PARAMETER 
370 MARK 4 *MARK THE TIME 
380 TRANSFER , TOEFL 
SS Oe 
400 * CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 

) 
411 GENERATE ae) Uae *CREATE THE FIRST 4 T-37'S 
412 ASSIGN onl “ASSIGN TO A PARAMETER 
413 TRANSFER DOWN 
414 GENERATE 120,10,60,4 CREAREMIINET IRS 4u—2 Ss 
415 ASSIGN 1,0 :ASSIGN TC A PARAMETER 
416 TRANSFER , DOWN 
417 * 
420 GENERATE 486,60, 1500 sCREATE T-37'S 
430 Weer NEGNRT.1,PISTA *THE SIMULATION STOPS AFTER 
AZ1 * "ALL THE A/C HAVE LANDED 
450 ASSIGN it -ASSIGN TO A PARAMETER 
440 TRANSFER , DOWN 
470 * 
480 GENERATE 480,60,1740 *CREATE T-2’'S 
490 TEST C NSGNRT,1,PISTA *THE SIMULATION STOPS AFTER 
491 * 7ALL THE A/C HAVE LANDED 
510 ASSIGN 1,0 :ASSIGN TO A PARAMETER 
520 DOWN #§ QUEUE CNTR *WAIT OUT OF THE RUNWAY 
530 GATE FV DANGER “DONT MOVE IF EMERGENCY HOLDS 
540 GATE SE CA *CHECK FOR A\C ON BASE LEG 
550 SEIZE CNTR -CAPTURE THE CONTROLLER 
560 DEPART CNTR *GOINS FOR LINE UP 

5 

S62 * LINE UP 
570 LINEUP ADVANCE TAs s LINE-UP CHECK 
580 RELEASE CNTR : TAKE-OFF 
596 MARK b *START FLIGHT TIMEFOR T-2 A/C 
600 MARI: 7 eSTART FLIGHT TIME FOR T-37 A/C 


af 





SCHED eee 


LE KR EF KE KHEKKEEK ERE KEKE KHER EE KEKE EEK KEE KEKHERKKEKKEKKKEE KER EKER KEKKKEEF FERRE 


20) 


% 


¥ 


= AIR-TRAFFIC CONTROL SIMULATION % 
BQ %& ee ee me ee ee ee % 
SQ x $e 
60 HHH HH HGH RH TE KRHKER ERK KKH KEK KHHE HEHEHE KHEKRHEKHEHHKHEKHKRRHHEAEKEKKHEHEHKHHERERE 
70 * 
Sa x 
70 x EXPONENTIAL DISTRIBUTION 
100 x 
110 EXP FUNCTION RNi C24 
O10/,.1 el Os meee a Sno ae 509/.5,-69/06,.9lo/.7,le2/e7oy 1s 38 
8,1.67.64,1,85/.88,2. 12/.9,2.5/.92,2. 92/94, 2.81/.95,2.99/ 70,502 
97 43.51.9855. -9/.99,4.6/.995 5.3/5 5d6 12 6.2/.999.7/.9998 28 
3 | a 
130 
140 MEAN FUNCTIGN 
Ly Vel)_ ele ete eos1s200. 488.6 
15006, 247. 63 7 23400, 557. 83/25200, 134° 5¢ 
218001494. 60/26600,221. 20 
150 *° 
140 # 
290 GENERST 9000 , -NEEXF sCREATE EMERGENCY EVENTS 
300 SEIZE DANGER 
310 ADVANCE 600, 180 sEMERGENCY HOLDS 
S20 RELEAS DANGER 
339 TERMINATE s EMERGENCY TERMINATES 
240 * 
350 GENERATE 9100, FNSEX® :4/C FROM OTHER BASES 
360 ASSIGN Le 3 ASSIGN TO A PARAMETER 
270 MARK -MARK THE TIME 
380 TRANSFER /TGEF 
390 % 
409 * CREATE AIRCRAFTS FROM TWO DIFFERENT SQUADROMS ON BASE 
10 # 
411 GENERATE 120,10, ,4 sCREATE THE FIRST 4 T-37°S 
412 ASSIGN i :ASSIGN TO A PARAMETER 
413 TRANSFER bo WN 
414 GENERATE 120, 10,60,4 sCREATE THE FIRST 4 T-2 
415 ASSIGN -ASSIGN TO A PARAMETER 
416 TRANSFER * DOWN 
%& 
420 GENERATE 600,60,1500 sCREATE T-37'S 
450 TEST N¢GNRT.1,PISTA >THE SIMULATION STOPS AFTER 
jase *ALL THE A/C HAVE LANDED 
450 ASSIGN 1,1 -ASSIGN TO A PARAMETER 
460 TRANSFER , DOWN 
AO x 
480 GENERATE 600,40, 1800 sCREATE T-2'S 
490 TESh N$GNRT.1,PISTA >THE SIMULATION STOPS AFTER 
491 * *ALL THE A/C HAVE LANDED 
510 ASSIGN 1,0 “ASSIGN TO A PARAMETER 
520 DOWN QUEUE CNTR *WAIT OUT OF THE RUNWAY 
530 GATE FV DANGER ?DONT MOVE IF EMERGENCY HOLDS 
540 GATE SE CAP6 *CHECK FOR A\C ON BASE LEG - 
550 SEIZE CNTR -CAPTURE THE CONTROLLER 
360 DEPART CNTR *GOING FOR LINE UP 
+ + 
562 * Woe see: 
570 LINEUP ADVANCE 70.20 sLINE-UP CHECK 
580 RELEASE CNTR : TAKE-OFF 
590 MARI z -START FLIGHT TIMEFOR T-2 A/C 
600 MARK 7 *START FLIGHT TIME FOR T-37 A/C 


08 
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TRAFFIC CONTROU SIMULATION 


KEKKEKKHKHK KEK KEEKEEK KE EKE KKH EK KHEK KEK KKH KH KEKKKEEK HHH KKH KH HK KH KHH HHH KHER EREHREE 


EXPONENTIAL DISTRIBUTION 


RN1 C24 

ERO Mich 6 0918/57 1.2/7.7 
RS ao eee Ae 20 Gt 69s 2.99/. 
9955. 3/.998,6.2/.999,7/.9998,8 


wy we 


‘a 
“ane 


De 
s ‘=’ 


ON 


— 
. 

—? 
_ 


C 
i> 


1,D8 
200,488.62 
=/25200, 124.58 


9000. FNSEXF *CREATE EMERGENCY EVENTS 
DANGER 

600, 180 s EMERGENCY HOLDS 

DANGER 


sEMERGENCY TERMINATES 


2A/C FROM OTHER BASES 
-ASSIGN TO A PARAMETER 
-MARK THE TIME 


9100, ENSEXE 
i? 


q 

, TOEF i 
FROM TWO DIFFERENT SQUADROMS ON BASE 
FNSMEAN , FNSEXF -CREATE AIRCRAFT 


oO AIR- 
40 +* ---- 
=O * 
60 

70 * 
80 * 
90 * 

100 # 

110 EXP FUNTTION 
ey. 104/.2,.222/.5,. 
eee 64,1.83/.88,2.12 
997s 5/.98,3.9/.99.4.6/ 
120'% 

130 # 

140 MEAN FUNCTION 
3600, 380. 11/5400, 221.60/ 
15000, 242. 63/23400,557.8 
31800, 494. 60/36600, 271. 20 
tao * 

160 +* 
290 GENERATE 

300 Seis 
310 ADVANC 
320 RELEASE 
330 TERMINATE 
340 * 
350 GENERATE 
360 ASSIGN 
370 MARK’ 
380 TRANSFER 
390 * 
400 * CREATE AIRCRAFTS 
AiG * 

420 GENERATE 

430 TEST L 

AzS TRANSFER 

SO ASSIGN 
460 TRANSFER 
470 +* 
S10 ACFTT2 ASSIGN 

S20 DOWN QUEUE 

30) Sie PY 

540 GATE SE 

550 SEIZE 

560 EPART 

570 LINUP ADVANCE 

S80 RELEASE 

590 MARK 

600 MARK 

610 ADVANCE 

615 TRANSFER 
620 +# 

630 SELECT E 
640 AREA TESTE 

650 +* 

660 * EAST AR 
670 +* 

680 SAVEVALUE 
690 SAVEVALUE 
700 WAIT QUEUE 


NSGNRT,1,PISTA >THE SIMULATION STOPS AFTER 
ee nee 

ef -ASSIGN TO A PARAMETER 

DOWN 

ale sASSIGN TO A PARAMETER 

CNTR *WAIT QUT OF THE RUNWAY 

DANGE -DONT MOVE IF EMERGENCY HOLDS 
CAP *CHECEK FOR A\C ON BASE LEG 

CNTR *CAPTURE THE CONTROLLER 

CNTR -GOING FOR LINE UP 

79,20 *LINE-UP CHECK 

CNTR > TAKE-OFF 

“START FLIGHT TIMEFOR T-2 A/C 

7 *START FLIGHT TIME FOR T-37 A/C 
100,10 SAFTER TAKE OFF TO DEPART POINT 
25, , OTHER -FORMATIONS, INSTR.FL., CPM 
2,101,108,0,F,QUEUP ;FIND EMPTY AREA IF EXISTS ONE 
BVSFIRST,1,WEST sCHOOSE EAST OR WEST AREA 
E Aes 

13+,P1 *RECORDS T-37'S 

14+,1 *RECORDS TOTAL # OF A/C 


sWAIT IN THE QUEUE 


-— 
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APPENDIX C 


This appendix contains the historical data for a three day 
period. For each aircraft the historical data consisted of the 
interarrival times for take off, the time spent waiting in the 
air traffic controller queue, the time spent waiting to ente: 


area 3 and the total flight time of T-37 aircraft. All times art 


stated in seconds. 


Interarrival times for take off 


INARRDAY! {Interarrival Times for take off Day 1) 

557 215 254 399 968 48 16 91 211 490 484 622 43 11 7 541 41 171 160 177 191 405 
472 62 27 646 205 1410 154 490 196 479 358 1637 207 B4 568 948 121) oll 
31 225 65 565 41 149 778 279 157 362 575 1147 305 288 475 SES 625 118 544 
274 1623 117 173 6 31 58 24 57 183 775 92 7 142 260 7 461 77 260 987 77 
29 378 225 235 766 70 183 200 79 322 22 99 118 298 70 3 280 146 


INARRDAY? (Interarrival Times for take off Day 2) 


291 226 259 597 89 70 351 1415 147 179 1006 368 167 43 478 22 135 18 21 287 207 
606 207 412 17 2Ad 959 414 182 1424 486 33 212 304 1168 362 261 477 2il 
$1 32 T10 14 128 267 304 90 37 274 11 B04 959 980 425 1491 280 887 6 230 


347 114 155 141 157 46 63 152 155 151 113 128 37 1101 1142 919 242 182 
1440 2L6 397 41 10829 392 26 221 BO 52 534 BS 99 261 S82 60 112 


INARSDAYS {Interarrival Tises for take off Day 3) 

432 170 1257 577 124 144 BB 378 260 634 ZIG 27 547 15 427 52 127 42 267 23 1k 
Q4 22% 144 596 118 798 421 148 171 657 214 20 371 106 1356 7B 129 22 al 
1.4 44 922 1135 481 1832 665 296 45 24 519 250 204 14 102 708 1258 1124 
1027 £76 564 914 221 49 33% 7 168 87 91 2 504 252 35 62 39 55 144 16d 67 
107 2 146 146 357 270 980 1422 869 334 1097 B73 277 447 312 621 274 2 al 
492 74 452 530 741 143 676 
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Time spent waiting in the air traffic controller queue 


Waiting Controller queue Area 3 
Time dayl day2 day3 dayl day2 day3 


O 42 45 49 af 
60 42 38 Sl 
120 12 20 11 


1 il 
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@ © See) @ Se) ei 
Oo0000rFc000OOR 
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COFOCOOFOOFOOFPW 


temeat flight time aircraft T-37 


Time dayl day2 day3 
3960 6 5 8 
4080 23 Ls 13 
4200 eZ 23 IL 
4320 10 9 8 
4440 4 3 2 
4560 0 Z a 
4800 0 al O 
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APPENDIX D 


This appendix contains the summary statistics and the 
hypothesis test results for the eight time periods for which 
the interarrival times for take off were homogeneous. These 


times periods were: 


1) 07230 = 705-40 5) 3405-14 oe 
Z) “O85420 —=20o710 6) 14:00 - 14:30 
3) 09710 sslls ao 7) 14:30 - 16:20 
A) SLi lO, eA 8). 26: 20G=. 17-230 


The Tables D.1 and D.2 contain the summary statistics 
for each time period and the figures D.1 through D.8 contain 
the distribution fitting and the Kolmogorov-Smirnov test 


results. 
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Table D.1 Summary Statistics 
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Ganple size 

AVerage 

Median 

Mode 

Geometric Mean 
VarLance 

Standard deviation 
Seannara error 

Man imcuon 
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Range 

Peat Guar lile 

Upper quartile 
aeerdiartile range 

a W162 55 

Standardized skewness 
Kurtosis 

Standardived kurtosis 
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i ee oe eee ee ane 
ee ORO RE! Onin e EE ¢ Oe EE RA OOD a Mr SMES OH SHS o 


sample S128 
Average 

Median 

100 2 

Geometric mean 
Var lance 
Standard deviation 
Standard error 
Manamum 

Ma imum 

Khange 

Lower quartile 

Upper quartile 
Interquartiie range 
Sk ewness 

Standardized skewness 
Purtosis 

Standardized kurtosis 
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INTERVAL 
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aie. 43 
=BO.114 oe \eaoe 7 
~60O 160 
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746.624 TS Fs pean 

120428 441355.7 
247.041 210.085 

58.6607 Shee 7G 

16 7 
1419 Ula e: 
eecy Ta 
144 42 
947 413 
4O% GVA 

1.650010 O.928297 
=.98766 2.04214 
Bays 78) —-O.118849 
229976 -O.162741 
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Figure D.1 Distribution Fitting Interval 1 
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Frequency histogram 
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Estimated EOLMOGOROV statistic DFPLUS = 0.124797 
Fsetinated KOLMOGOROV statistic DMINUS = 0.111456 
Estimated overall statistic DN = 0.124797 


Approximate significance level = 0.286891 


Figure D.2 Distribution Fitting Interval 2 
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Approximate significance level 


Figure D.3 Distribution Fitting Interval 3 
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Fetimated KOLMOGOROV statistic DFLUS = 6,0935019 
Estimated FOLMOGOROYV statistic DMINUS = 0,.0871474 
Estimated overall statistic DN = 0.09339019 
Approximate saL.ynifrcance level = 1 


Figure D.4 Distribution Fitting Interval 4 
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Frequency Histograti 
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ERsetimated EFOLMOQGORDV statistic DFLUS = O.0625072 
Estimated FOLMOQGOROV statistic DMINUS = @.,110322 
Estimated overall statistic DN = O,110222 
Approximate Significance level = 0.999941 


Figure D.5 Distribution Fitting Interval 5 
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Frequency Histogran 
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NTERVALG 


Estimated FOLMOGOROV statistic DFLUS = 0.128466 
Estimated FOLMOGB0ROV statistic DMINUS = O.0882622 
Estimated overall statistic DN = 0.128466 
Approximate significance level = 9.940454 


Figure D.6 Distribution Fitting Interval 6 
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Estimated FOLMOQGOROV statistic DFLUS = O.,O35223982 
Retimated KOLMOGOROV statistic DMINUS = 0.117 
Estimated averall statistic DN = G.117 
Approximate siqynificance level = 0.999736 


Figure D.7 Distribution Fitting Interval 7 
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Fstimated KOLMOBORDY statistic DFLUS = 0,90867169 
cetimated KOLMOGOROV statistic DMINUS = GO,0760762 
Estimated overall statistic DN = 0.08467169 
Approximate significance level = 1 


Figure D.8 Distribution Fitting Interval 8 
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